Feshbach Resonance Cooling of Trapped Atom Pairs 
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Spectroscopic studies of few-body systems at ultracold temperatures provide valuable information 
that often cannot be extracted in a hot environment. Considering a pair of atoms, we propose a 
cooling mechanism that makes use of a scattering Feshbach resonance. Application of a series of 
time-dependent magnetic field ramps results in the situation in which either zero, one, or two atoms 
remain trapped. If two atoms remain in the trap after the field ramps are completed, then they 
have been cooled. Application of the proposed cooling mechanism to optical traps or lattices is 
considered. 

PACS numbers: 32.80.Pj, 03.75.-b, 34.50.-s 
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A Feshbach resonance |H 13- El occurs for two atoms 
when their collision energy becomes degenerate with a 
vbound state in a closed collision channel, producing brief 
transitions into and out of this state. In recent years, 
these resonances have been used extensively to control 
the interaction strength in dilute atomic gases 0, [H, El • 
Here, we utilize some of the unique characteristics of a 
Feshbach resonance to develop a cooling mechanism that 
is applicable to two externally-conhned atoms. Applica- 
tion of a series of field ramps (i.e., cooling cycles) leads to 
a cooled atom pair (an atom pair with reduced internal 
energy), provided that the atom pair remains trapped. 

In the following, we first develop the basic mechanism 
of the Feshbach resonance cooling process. The feasi- 
bility and effectiveness of the proposed scheme are then 
illustrated through an application to a realistic system 
of two atoms in a trap. Finally, possible applications to 
optical traps are discussed. 

The concept of Feshbach resonance cooling grows out 
of the observation that the quantum-mechanical energy 
levels of two atoms in a harmonic trap shift by an en- 
ergy corresponding to approximately two trap quanta, 
as a control parameter is swept in one direction across 
the resonance. Throughout this article, we refer to this 
control parameter as the magnetic field B used to manip- 
ulate the atom-atom scattering length a in the vicinity 
of a pole. In other contexts, the shift of the energy lev- 
els could be introduced by varying the detuning of an 
off-resonant dressing laser, or by varying an electric field 
strength. The ideas presented here in terms of the con- 
trol parameter B can be straightforwardly extended to 
those other contexts. 

The Schrodinger equation for two interacting identi- 
cal mass m atoms under spherical harmonic confinement 
with trapping frequency v decouples into two equations: 
one involving the three relative coordinates of the pair, 
and another involving the three center-of-mass (CM) co- 
ordinates 0, • We consider the Schrodinger equation in 
the relative coordinate for two trapped atoms interacting 
through a central potential and assume for the time be- 
ing that the center of mass coordinate is translationally 



cold. Accounting for an applied external magnetic field 
B through a independent quantum defect (3ei(B), the 
energies E n i(B) associated with the relative motion of 
an atom pair are given by Q 

E nl (B) = (2n-2p El (B)+l + 3/2)hLJ, (1) 

where lo = 2txv. Here, the quantum defect Pei(B) de- 
pends strongly on the relative orbital angular momen- 
tum I of the pair, while it depends only weakly on the 
radial oscillator quantum number n. The dependence of 
Pei(B) on the energy is weak on the scale of an oscillator 
quantum, i.e., \d/3Ei(B)/dE n i\ <C l/fiw. 

As will become clear later, the quantum defect for 
one relative partial wave I for an atom pair, e.g., the 
s-wave, p-wave, or cZ-wave, must rise by unity across the 
energy range k B AT of interest, and across the accessi- 
ble range of the control parameter, AB. In fact, this 
variation of Pei(B) by unity corresponds to the Fesh- 
bach resonance, which causes the scattering phaseshift 
to rise by tt. A simple closed-form expression exists 
for (3 ei H,IS, which simplifies at energies higher than a 

few trap quanta to @ei(B) ~ arctan ( %~j^jf~~^ where 
a(E n i,B) is the energy- and field-dependent scattering 
length and L osc — ^Jh/ (/mj) with /i = m/2 denotes the 
characteristic oscillator length. 

In this paper, we focus on an s-wave resonance, though 
this formalism can be readily extended to higher partial 
wave resonances. When an s-wave Feshbach resonance 
occurs, the limiting low-energy scattering phaseshift 
is proportional to the wavenumber k — (2fj,E/h 2 ) 1 ^ 2 . 
Omitting the subscript I, the E- and B-dependent scat- 
tering length is then given by 
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where at, g is the background scattering length. At the 
magnetic field strength B rcs of the resonance a zero- 
energy bound state occurs. The resonance width in en- 
ergy is related to the width in the control parameter 
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FIG. 1: (Color online) Energy levels E n for the relative co- 
ordinate of a harmonically trapped 85 Rb atom pair near the 
Bo ~ 155.2 G Feshbach resonance, as a function of the mag- 
netic field B. A rather large trapping frequency of v = 1 MHz 
is used in order to clarify the field dependence of the energy 
levels. Cooling is performed by ramping the magnetic field B 
slowly from B\ to B2 and then quickly back to Si. A more 
realistic ramp would likely encompass more level curves (i.e., 
cover a larger field range) . The state which undergoes a shift 
for B = B2 (which we will label as n — Q later) is indicated 
by a dashed line. Ideal cooling is described diagramatically in 
the inset (same axes), where population transfer from point 
a to point b occurs during the slow field ramp and from b to 
c during the fast ramp. See the text for details. 



A by T E — 2kab g E; es (B)A, where E' res denotes the rate 
at which the resonance energy E ICS varies with the con- 
trol parameter [Tjj • Figure[Jillustrates the characteristic 
s-wave energy levels E n appropriate for the relative mo- 
tion of two atoms in a spherical harmonic oscillator trap, 
as functions of the applied magnetic field B for a mag- 
netic Feshbach resonance in 85 Rb(2, -2)+ 85 Rb(2, -2). 
In Fig. ^ the parameters adopted are B YCS = 155.2 G, 
E' ICS = -3.5MHz/G, T B = 10 G, and a bg = -380 a , 
where ao is the Bohr radius. 

Feshbach resonance cooling entails ramping the mag- 
netic field through the region where the energy levels 
shift by w 2fhu. Figure ^ denotes the initial B-field by 
B\ . For an atom pair taken from a source of atoms with 
temperature T, the probability of the pair being in the 
energy eigenstate \n(Bi)), and hence, of having the en- 
ergy E n (Bi), is given by the corresponding Boltzmann 
factor. By ramping the magnetic field from B\ to B2 
(also shown in the inset of Fig.QJ, sufficiently slowly to 
be adiabatic, the population of each energy level remains 
unchanged, while the energy level itself is reduced com- 
pared to its value at B-y, E n (B 2 ) ~ E n (Bi) — 2hui; the 
atom pair has lost energy. To further reduce the energy 
of the two atoms, another field ramp has to be applied. 
In order to do this, it is necessary to return the mag- 
netic field to B\ without adding the energy back that 
has been just removed. A fast, nonadiabatic change of 
the magnetic field from B 2 to Bi ensures this, since this 
change simply projects the eigenstate \n(B2)) onto the 



eigenstate \n(Bi)). Note that this proposed energy re- 
duction (cooling) comes at the expense of an energy gain 
of a single level (or a few levels, see below for a detailed 
discussion). 

To model the effects of the magnetic field ramps, we 
have developed a two-channel Feshbach resonance model, 
based on the single-channel model described in Ref. ^lj • 
Both of these models describe a two-atom Feshbach res- 
onance for a harmonic trap, and utilize a zero-range po- 
tential to describe the interaction between the two atoms. 
The two-channel model has the advantage of allowing for 
a field-dependent resonance state. In the two channels, 
the s-wave radial solutions for the relative coordinate r 
of the atom pair satisfy the equations 
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-^-^ + ^r 2 ju 2 {r) = (E-£)u 2 (r), (4) 

where e is the energy shift of the second channel from 
the first channel. The zero-range potential imposes a 
boundary condition at the origin, which is parameterized 
as 
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A quantum-defect-theory treatment, similar to Ref. |llj |. 
can then be applied. The scattering length predicted by 
this model (when to — > 0) is 
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a(E,B)= - 
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which can be compared to the measured scattering length 
to determine the values of the parameters a±, a 2l and /?. 
The parameters also affect the magnetic-field dependence 
of the adiabatic energy states, and their adjustment is 
able to provide satisfactory agreement with experimen- 
tal data in the regions of interest to us. E.g., for 85 Rb, we 
find ai = -435 a , a 2 = 1.485 a , and = 0.0011618 a^ 1 . 
Simulations can then be performed by specifying an ini- 
tial state of the system and numerically propagating the 
Schrodinger equation. 

The simulations reveal the effect of the adiabatic and 
nonadiabatic field ramps: Assume first that the atom 
pair is in a pure state at B = B\. As expected, the adia- 
batic field ramp (B\ to B 2 in Fig. decreases the energy 
of the atom pair irrespective of the initial eigenstate cho- 
sen. A nonadiabatic ramp (B 2 to Bi in Fig. causes 
a state at B 2 that is not degenerate with the resonance 
state to project onto a state at B±, with approximately 
the same energy as the initial state at B 2 . However, if 
the pair is initially in the state at B 2 that is degener- 
ate with the resonance state, the fast ramp results in a 
strong projection onto the resonance state. In this case, 
the atom pair gains energy since the resonance state at 
Bi has a higher energy than the initial state at B 2 . 



We now generalize our scheme to the more practical 
situation of a mixed initial state, and show how removal 
of hot atoms leads to an intriguing cooling scheme. For a 
mixed state, the occupation probability of a pair level 
with energy E n in the relative motion is determined 
in terms of a Boltzmann factor by er En l T /Z(t) with 
r = k B T, where k B is Boltzmann's constant, Z(t) is the 
(relative) partition function, and T is the temperature 
of the source of the two atoms. Based on the results of 
field ramps for pure states discussed above, we see that 
application of a cooling cycle (slow ramp from Bi to B 2 
plus fast ramp back to B\) for a mixed state will do two 
things: (a) decrease by 2ftw the energy of the population 
in states which undergo a full energy shift between B\ 
and B2, and (b) increase the energy of the population in 
the state that is degenerate with the resonance (i.e., un- 
dergoing an energy shift) at B2 by moving it to the state 
(or states) degenerate with the resonance at B\ . We will 
denote this state from which the heated population orig- 
inates as n = Q. 

Figure |21 illustrates the effect of a single cooling cy- 
cle on a mixed state, using the same Feshbach resonance 
as shown in Fig. ^ The black line represents an ini- 
tial s-wave probability distribution for the states associ- 
ated with the relative coordinate of an atom pair in a 
harmonic trap with v = 1 MHz and source temperature 
T = 1 mK. The red line represents the same probability 
distribution after application of a slow and a fast mag- 
netic field ramp, for Q = 10. Application of one cooling 
cycle moves the population of the state Q to states with 
much higher energy, here n « 85, evidenced by the spike 
in Fig. |3 At the same time, the field ramps move the 
population of each state with n > Q to the next-lowest 
state, which has ~ 2hu> less energy. Our numerical sim- 
ulations indicate that, on average, the net energy of the 
system is increased after application of one cooling cycle. 
The few cases where the pair gains a large amount of en- 
ergy overcomes the many cases where the pair looses a 
small amount of energy. This behavior is expected, since 
Ketterle and Pritchard 0] m 1992 pointed out the im- 
possibility of creating a cooling scheme relying solely on 
time-dependent potentials. 

An estimate of the efficiency of cooling can be made by 
assuming that the population of the state Q is removed 
from the trap, while the population of all states with 
n > Q are moved to the next-lowest state, that is, to 
states with n — 1. This assumes that the range of the field 
ramps is such that the heated fraction (n = Q) ends up at 
an energy corresponding to negligible thermal population 
(this is the case in Fig. [2J , and that all population above 
a specified energy can be removed. If we approximate 
the level energies at B = B\ by E n {B\) w 2nHu>, the 
probability to remove an atom pair during a cooling cycle 
is 



-2HloQ/t 
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FIG. 2: (Color online) Illustration of the effect of a single Fes- 
hbach resonance cooling cycle for T = 1 mK and v = 1 MHz. 
The black (red) line represents the population distribution 
before (after) application of one slow and one fast magnetic 
field ramp. The state Q (here Q — 10) is indicated. Popula- 
tion from n — Q and nearby states is moved to higher states 
with n ~ 85. Population initially in a state with n > Q is 
moved to the next-lowest state (see the inset close-up). 



the energy of the n ~ Q population removed from the 
trap, plus the energy loss for states n > Q: 



AE(Q,t) = (2huQ + (E CM )) 
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Noting that {Ecu) — 3t (since (E to t) — 3r for a single 
atom in a harmonic trap), and with X)^Lq+i e~ 2hu)n l T m 



-2huQ/r 



r/2hu>, Eq. (JSJ becomes 

-2HujQ/t 

AE(Q, r) = - {2hwQ + At) . 
Z(t) 



(9) 



The energy efficiency E c g, defined as the amount of en- 
ergy removed per atom removed, is then given by 



£cff(<2, r) = 2hioQ + At. 



(10) 



Since Q determines the efficiency of the cooling process, 
it is referred to as the cooling parameter. Results from 
our numerical model indicate that Eq. (|10f> provides a 
good estimate of the efficiency. 

The time scale for one cooling cycle is determined by 
the speed of the adiabatic field ramp. This speed in turn 
is determined by the strength of the coupling between the 
resonance state and the trap states. The smaller the cou- 
pling for an avoided crossing, the slower is the field ramp 
required to maintain adiabaticity. The coupling between 
the resonance state and the trap states is related to the 
resonance width parameter Te, which can be used in a 
Landau-Zener estimate of the transition probability |lOj | , 



The average energy decrease in a cooling cycle is due to 



P tr = exp - 



\dB/dt\ \dE/dB\) 



(11) 
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FIG. 3: Probability that a pair of atoms remains trapped 
vs. the average total kinetic energy of the two atoms in os- 
cillator units (note that k B T = {Eu>t)/6 for two harmoni- 
cally trapped atoms). Three different cooling parameters are 
used: 2huiQ = 5r (solid line), 9r (dashed line), and 12r (dot- 
dashed line). It is assumed that rethermalization occurs be- 
tween cooling cycles (see text). Inset: probability to remain 
trapped vs. the number of cooling cycles for the same three 
cooling parameters. 

Motivated by the possibility of experimentally trap- 
ping a small, deterministic number of atoms |13| . we 
now explore the experimental feasibility of our cooling 
scheme. A Feshbach resonance cooling experiment in- 
volves a sequence of cooling cycles. As discussed above, 
a single experiment could result in a heated atom pair, 
which in turn would be lost from the trap. To end up with 
a cooled atom pair, multiple ramp cycles are required. To 
see the effect of cooling, Eqs. (J7J and (|SJ) can be iterated. 
For a variety of cooling efficiency parameters Q, Fig. 
shows the probability for an atom pair to remain trapped 
vs. the average total kinetic energy (the energy of both 
the relative and the CM degrees of freedom) of the two 
atoms in oscillator units. Included in this calculation is 
the probability that the atom pair is in an s-wave state 
to begin with, because the field ramp has no effect on 
other partial waves. We assume that rethermalization 
occurs between cooling cycles, which could be ensured 
by, for example, introducing a slight anharmonicity into 
the trapping potential. 

To be more specific, we consider a crossed-beam op- 



tical dipole trap [T^ which offers a good blend of large 
trap frequency (for a large s-wave fraction) , isotropy, and 
anharmonicity (for rethermalization between the relative 
and CM degrees of freedom). Assuming the dipole trap 
has an average frequency of v = 10 kHz and contains two 
atoms taken from a source with temperature T = 8 /jK 
{{Etot) / Tkjj = 100), we see from Fig. [3J solid line, that a 
temperature of 0.16 /iK ((E to t)/fujo = 2, both atoms in 
the ground state) could be reached 10% of the time by 
performing less than 20,000 cooling cycles. For a range 
in magnetic field for the ramps of AB ks 1 G, and us- 
ing Eq. Ijllfl with Pt r = 0.1, we see that such a series of 
field ramps could take place in under 1 s. A pertubative 
calculation accounting for the trap anharmonicity (see 
Ref. for details of a similar treatment) indicates that 
rethermalization between the relative and CM degrees of 
freedom should occur on a time scale comparable to a sin- 
gle ramp time for a crossed-beam dipole trap. This will 
ensure that the relative s-wave distribution will rether- 
malize with each ramp and that the cooling of the rela- 
tive coordinate will also cool the CM coordinate (both of 
which we have assumed up to this point). 

Another possibility would be to apply our cooling 
scheme to atom pairs trapped in an optical lattice. In 
this case, field ramps could be performed on the lattice 
ensemble of atom pairs, with a certain percentage of sites 
resulting in cooled pairs, while other sites will have either 
zero atoms or one (uncooled) atom. It may also be possi- 
ble to prepare the optical lattice by some other means to 
have a high probability of exactly double occupancy at 
each lattice cite (see for example Ref. |r| ) • From such an 
initial state, a Feshbach resonance cooling scheme could 
be used to efficiently cool atom pairs to low-lying trap 
states. 

In summary, we have proposed a novel cooling scheme 
that makes use of Feshbach resonances. We have also 
shown that such a cooling scheme may be experimen- 
tally feasible, although it may not be competitive with 
proven cooling methods like evaporation. At the same 
time, plenty of room for exploration remains. Extension 
to atom clouds might also be possible. 
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